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利用 贝 叶 斯 分 层 模 型 估算 炊 变 体 的 多 普 勒 因子 


宁 宫 明 ， 黄 伟 天 ， 毛 李 胜 、 


(云南 师范 大 学 物理 与 电子 信息 学 院 ， 云 南 昆明 650500) 


摘要 : 基于 亚 毫 米 波 阵列 (Submillimeter Array, SMA ) 的 1 mm 波长 的 长 时 间 监 测 数据 ， 
利用 贝 叶 斯 分 层 模 型 对 沭 变 体 的 光 变 曲线 进行 拟 合 ， 估 算 了 155 PH R AMA BE Fe 
变 多 普 勒 因子 。 利 用 Wilcoxon 秩 和 检验 ， 比 较 了 蝎 虎 座 BL 型 天 体 (BL Lac) 子 样本 和 平 谱 射 
电 类 星体 (Flat Spectrum Radio Quasar, FSRQ) 子 样本 的 亮 温度 和 多 普 勒 因子 分 布 ， 也 比较 了 
费 米 焰 变 体 和 非 费 米 浪 变 体 的 亮 温度 和 多 普 勒 因子 分 布 。 结 果 表 明 ;: (1) 相对 蝎 虎 座 BL 型 
天 体 ， 平 谱 射 电 类 星体 平均 具有 更 高 的 多 普 勒 因子 ; (2) 相 对 非 费 米 焰 变 体 ， 费 米 漆 交 人体 平 
均 具有 更 高 的 亮 温 度 和 多 普 勒 因子 。 通 过 与 15 GHz 波段 的 亮 温 度 和 多 普 勒 因子 进行 比较 发 
现 ， 射 电 亮 温度 和 多 普 勒 因子 可 能 有 向 更 高 频率 降低 的 趋势 。 
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活动 星系 核 ( Active Galactic Nucleus, AGN ) 是 宇宙 中 电磁 辐射 最 强 的 天 体 之 一 ， 一 部 分 活动 星系 
核能 产生 向 外 运动 的 喷 流 并 发 出 同步 辐射 '" 。 当 活动 星系 核 的 喷 流 方向 与 观测 视线 非常 接近 时 ， 这 
类 源 称 为 次 变 体 ( Blazar) ， 其 主要 观测 特征 有 高 光度 、 快 速 光 变 、 高 偏振 、 视 超 光 速 运动 等 。 为 了 
解释 这 些 特殊 现象 ， 文 [2] 提 出 了 一 个 活动 星系 核 的 统一 模型 。 在 该 模型 中 ， 炊 变 体 中 心 有 一 个 超大 
质量 黑洞 ， 周 围 围绕 着 一 个 吸 积 盘 ， 喷 流 从 黑洞 中 心 以 垂直 于 吸 积 盘 的 方向 喷 出 。 兆 变 体 的 多 波段 辐 
射 均 由 于 多 普 勒 效应 而 被 放大 ' 沾 。 根 据 光 谱 中 的 发 射线 性 质 ， 炊 变 体 可 以 分 为 蝎 虎 座 BL 型 天 体 和 平 
谱 射 电 类 星体 两 个 子 类 。 前 者 光谱 中 没有 发 射线 或 仅 有 弱 发 射线 ， 后 者 表现 为 强 的 发 射线 六 。 

量化 喷 流 的 参数 之 一 是 多 普 勒 因子 6 = [LT(1 -Beos0] =l I 其 中 万 是 洛 伦 兹 (Lorentz) 因子 ， ELH 
三 =1IZV(CL-B)，B 是 以 光速 为 单位 的 喷 流 运动 速度 ，6 是 喷 流 方向 与 观测 视线 的 夹 角 '5。 但 是 ， 
无 论 是 6 还 是 9 都 难以 通过 测量 直接 得 出 ， 所 以 直接 确定 多 普 勒 因子 较为 困难 。 考 虑 到 多 普 勒 因子 对 
研究 炊 变 体内 豪 辐 射 的 重要 性 和 必要 性 ， 研 究 人 员 提 出 了 许多 间接 估算 多 普 勒 因子 的 方法 。 文 [6] 假 
设 观测 到 的 X 射线 来 自 同步 辐射 光子 的 逆 康 普 顿 散射 过 程 ， 通 过 比较 观测 X 射线 流量 和 根据 其 长 基 
线 干涉 测量 ( Very Long Baseline Interferometry, VLBI) 观测 数据 的 预言 X 射线 流量 ， 定 义 道 康 普 顿 多 普 
勒 因 子 。 文 17] 假 设 噶 流 中 磁场 和 辐射 粒子 能 量 均 分 ， 比 较 均 分 亮 温 度 和 观测 亮 温 度 ， 从 而 得 出 源 的 
均 分 多 普 勒 因子 。 文 [8] 和文 [1] 通 过 将 射电 光 变 曲线 分 解 为 指数 上 升 和 指数 下 降 爆 发 过 程 ， 确 定 爆 
发 亮 温 度 ， 最 后 对 比 均 分 亮 温 度 来 估算 焰 变 体 的 多 普 勒 因子 ， 称 之 为 光 变 多 普 勒 因子 。 文 [9] 利 用 贝 
叶 斯 分 层 模 型 ， 将 光 变 曲线 建 模 为 随机 背景 上 的 指数 上 升 和 指数 下 降 过 程 ， 通 过 与 均 分 亮 温 度 比 较 估 
算 了 837 个 泡 变 体 的 光 变 多 普 勒 因子 。 文 [10] 通 过 拟 合 炊 变 体 的 多 波段 能 谱 分 布 (Spectral Energy 
Distribution，SED ) 估 计 了 999 个 炊 变 体 的 多 普 勒 因子 。 文 13] 提 出 了 一 种 基于 噶 流 辐射 功率 和 宽 线 区 
光度 之 间 的 相关 性 ， 佑 算 了 350 个 光 变 体 的 多 普 勒 因子 。 以 上 几 种 多 普 勒 因子 估算 方法 均 基 于 不 同 的 
假设 ,对 于 单个 源 而 言 ， 很 难 确定 方法 之 间 的 优 劣 或 准确 性 。 
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本 文 利用 亚 毫 米 波 阵 列 1 mm 波长 的 长 时 间 监 测 数据 ， 采 用 与 文 [ 9] 相同 的 方法 ， 佑 算 了 155 ARE 
变 体 的 射电 亮 温 度 和 多 普 勒 因子 ， 则 在 比较 次 变 体 子 类 的 亮 温 度 和 多 普 勒 因子 分 布 差异 。 本 文采 用 的 宇 
宙 学 参数 ， 甩 =71 km ' s + Mpe, 2,,=0.27° 1 。 无 特殊 说 明 ， 文 中 亮 温 度 取 以 10 为 底 的 对 数值 。 


1 样本 和 数据 


亚 毫米 波 阵列 位 于 夏威夷 莫 纳 克 亚 山峰 。 它 的 一 项 重要 科学 任务 是 持续 监测 河 外 射电 源 :21 ， 监 
测 波长 为 mm 和 850 hm， 前 者 的 监测 密度 远大 于 后 者 。 亚 毫米 波 阵列 的 流量 数据 定期 更 新 ， 授 权 使 
用 后 可 通过 SMA 官网 下 载 (http ://smal.sma.hawaii.edu/callist/callist.html) 。 通 过 与 第 5 版 罗 CI ANN 
星 表 ( Roma-BZCAT) 进行 小 于 2" 的 位 置 匹配 52] ， 亚 毫米 波 阵 列 目 前 共 监 测 了 436 个 耀 变 体 。 通 过 肉眼 
仔细 检视 这 些 源 的 光 变 曲线 ， 我 们 挑选 155 个 在 1 mm 波段 具有 明显 爆发 现象 旦 采样 相对 密集 的 次 变 
体 。 这 个 样本 包括 25 个 蝎 虎 座 BL 型 天 体 ，112 个 平 谱 射电 类 星体 ，18 个 类 型 未 知 的 耀 变 体 (BCUs ) 。 
另外 ,通过 与 第 4 版 费 米 活动 星系 核 星 表 (4LAC) 匹 配 '' ， 我 们 发 现 样本 中 的 104 个 耀 变 体 具有 伽 马 
射线 辐射 ， 称 为 费 米 炊 变 体 ( Fermi blazars) ， 其 余 的 源 称 为 非 费 米 炊 变 体 (Non-Fermi blazars ) 。 作 为 示 
例 ， 图 1 展示 了 样本 中 两 个 次 变 体 的 1 mm 光 变 曲线 。 
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图 1 AAHH 1 mm 光 变 曲线 


Fig. 1 1 mm light curves of two blazars 


2 J 法 

我 们 采用 与 文 [9] 相 同 的 方法 对 155 SERA 1 mm 光 变 曲线 进行 拟 合 。 本 文 仅 简 要 描述 该 方 
法 ， 具 体 方法 细节 和 Python 程序 代码 可 见 文 [9] 、 文 [15] 和 https ://github.com/dhuppenkothen/magnetron2/ 
tree/blazars。Magnetron 是 一 个 用 Python 语言 实现 的 贝 叶 斯 分 层 模 型 ， 它 将 光 变 曲线 建 模 为 一 系列 爆 
发 过 程 的 晋 加 ， 每 个 爆发 过 程 分 为 指数 上 升 过 程 和 指数 下 降 过 程 。 爆 发 过 程 的 个 数 是 一 个 自由 参数 ， 
每 个 爆发 过 程 由 4 个 参数 表征 ， 即 位 置 、 光 变 幅度 (以 Jy 为 单位 ) 、 上 升 时 间 ( 以 天 为 单位 ) 和 偏 度 
( 即 下 降 时 间 与 上 升 时 间 之 比 ) 。 光 变 幅 度 定义 为 峰值 流量 和 背景 水 平 之 间 的 差 值 。 图 2 显示 了 利用 
Magnetron 贝 叶 斯 分 层 模 型 拟 合 两 个 源 1 mm 光 变 曲线 的 结果 。 

利用 文 [9] 的 方法， 我 们 采用 公式 


LAS, (r) 
vt, (1 +z) 
估算 亮 温度 ， 其 中 ，z 为 红 移 ; ASw(z) ADERIR MN Jy); di 为 光度 距离 (单位 为 Mpc); vH 
观测 频率 (单位 为 CHz) ; ,为 爆发 的 上 升 时 间 ( 单 位 为 天 )。 人 针对 每 个 爆发 过 程 计算 7,,,， 并 找到 其 中 


T. = 1.47 x 101 (1) 
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的 最 大 值 7,,,,。。 重 复 上 述 过 程 ， 对 所 有 可 用 的 模型 (平均 有 157 个 ) 创建 7,, ,的 分 布 。 亮 温度 的 计 
算 结 果 及 误差 取 7,, ,分 布 的 中 值 和 1o 误差 。155 个 泡 变 体 中 ， 有 两 个 蝎 虎 座 BL 型 天 体 没有 红 移 ， 
在 计算 亮 温度 时 取样 本 中 蝎 虎 座 BL 型 天 体 的 红 移 平均 值 。 图 3 展示 了 两 个 源 的 了 ,分布 图 。 
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图 2 两 个 源 的 观测 ( 黑 点 ) 和 后 验 (红线 ) 取 样 光 变 曲线 。 蓝 色 虚 线 表示 分 解 出 的 独立 爆发 过 程 
Fig.2 Observed (black points) and posterior sampled (red lines) light curves of two sources. 


The blue dashed lines represent the decomposed independent flares 
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图 3 两 个 源 的 亮 温 度 了 分 布 ， 红 色 虚 线 表 示 亮 温度 分 布 的 中 位 数 ， 蓝 色 和 绿色 之 间 的 部 分 表示 lo 置信 区 间 


Fig.3 The brightness temperature T» ma distribution of two sources, where the red dashed line represents the median brightness 


temperature and the region between the blue and green dashed lines represents the lo confidence interval 


泡 变 体 的 辐射 特性 与 聚 束 效应 有 关 ， 应 考虑 多 普 勒 因子 修正 。 我 们 利用 文 [9] 的 方法 计算 了 光 变 


多 普 勒 因子 ， ' 
T... 
Ô ar =(1 +z) T + (2) 


为 了 计算 5,,,， 我 们 从 每 个 源 的 7,, ,分布 的 中 值 和 Lo 误差 范围 内 取出 一 个 随机 值 。 根 据 文 [9] ， 从 
平均 值 为 (7.,〉》= 2. 78 x 10"K， 标 准 偏差 为 a， = 7.2 x 10" 的 高 斯 分 布 中 取出 7 的 随机 值 。 我 们 使 
用 (2) 式 计算 光 变 多 普 勒 因子 8,.。 重 复 上 述 过 程 1 000 次 ， 从 而 为 每 个 源 创建 了 6,, 的 分 布 。 最 后 ， 
根据 每 个 源 的 6, 分布 ,得 到 多 普 勒 因子 的 中 位 数 和 lc 置信 区 间 。 


3 ”结果 与 讨论 
本 文 共 估算 了 155 个 耀 变 体 在 1 mm 波段 的 亮 温度 和 光 变 多 普 勒 因子 ， 结 果 见 表 1。 第 1 列 为 源 
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名 ; 第 2 列 为 次 变 体 类 型 ， 第 3 列表 示 源 是 否 为 费 米 炮 变 体 ，Y HAKR, N WIERA REEE; 
第 4 列 为 红 移 ; 第 5 列 为 对 数 亮 温度 及 其 误差 ,单位 为 K; 第 6 列 为 光 变 多 普 勒 因子 及 误差 。 图 4 和 
图 5 展示 了 155 个 炊 变 体 的 射电 亮 温 度 和 光 变 多 普 勒 因子 的 分 布 情况 。 
表 1 155 个 耀 变 体 的 射电 亮 温 度 和 光 变 多 普 勒 因子 
Table 1 The radio brightness temperatures and variability Doppler factors of 155 blazars 


Name Type Fermi z logT,,,/K Svar Name Type Fermi z logT,,,/K Svar 
(1) (2) (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 
J0006-063 BL Lac N 0.347 12153 3.25%) | J0457+067 FSRQ 


= 


0.405 11.422 2 3721 


J0019+734 FSRQ 1.781 13.3417 19. 76°28? | J0457-234 FSRQ 1.003 13.35220 12. 71°78 


J0050-094 BL Lae 635 10.74172 1.785 |J0501-019 FSRQ 2.291 13.66 34.42225 


J0102+584 BCU 644 13.0913 9.773 | JO510+180 FSRQ 0.416 13.215 7.145% 


J0106-405 FSRQ 584 12.8222% 10.01% || Jo522-364 BCU 0.056 12.8453 4. 827703 


J0121+118 FSRQ 57 11.3677 2.62 |] JO530+135 FSRQ 2.07 13.91090 32.311513 


J0136+478 FSRQ 859 13.77209 14.61157 || J0532+075 FSRQ 1.254 13.12 13.56123 


° ° e e ° = 


J0137-245 FSRQ .835 12.73214 9.73: ||J0533+483 FSRQ 1.16 11.8375 3.79 


= 


J0204-170 FSRQ .74 12.23209 8.77:%3 | J0538-440 BL Lac .892 13.87229 16.4837 


J0205+322 FSRQ 


= 


.466 12.73232 9.88: || J0607-085 FSRQ 


-1.18 


0 
0.87 108208 2.6470 
0 


J0224+069 FSRQ .511 14. 15712 18.09t$2 || J0609-157 FSRQ .324 13.2178 7.051 


o Oo 


J0228+673 FSRQ .523 11.367 2% 2. 46°55, || J0646+448 FSRQ 3.396 10. 14°73 4. 46°53 


J0237+288 FSRQ 


= 


.206 13.01725 12. 03$; 5 || Jo725+144 FSRQ 1.038 14. 58°13) 42.7954% 


J0238+166 BL Lac 94 13.87106 15. 12°38 | J0730-116 FSRQ 1.589 12.5254 7.712 


J0241-082 BCU .005 8.587792 0.2273 |J0738+177 BL Lac .424 10.87777 2.35778 


No ° 


68 107427 3.8" || JO739+016 FSRQ .189 13.37 6.962 


J0303+472 BL Lac 11.59%; 2.421 || JO748+240 FSRQ 


= J0309+104 FSRQ 863 12.92: 9.243% |J0750+125 FSRQ 


0 
0 
0.409 13.15%], 9.115 
0.889 12.4371 6. 92°39 
0 


JO310+382 FSRQ 945 13.02113 15.49°7% || J0757+099 BL Lac .266 14.11 15. 2874 2 


JO313+413 BCU 136 12.42227 4.931172 || JO808-078 FSRQ 1.837 14. 017353 33. 01731, 


° ° ° ° 


J0319+415 BCU .018 9.5157 0.312353 | JO825+031 BL Lac 506 13.79114 13.46552 


J0334-401 BL Lac 


= 


.445 13.11255 17.13127 | J0826-225 BL Lac 122 5.546 


= 


J0336+323 FSRQ .259 13.5811% 15.89172 | J0830+241 FSRQ 939 14.32120 23. 8777 


Pb ° ° = 
Kej 
= 


J0339-017 FSRQ .85 13.59: 14. 36°38 | J0841+708 FSRQ .218 14.29 37.46156 


(= =. 


J0348-278 FSRQ .991 13.1943% 14.7712% | J0854+201 BL Lac 


i>) 
Ww 
S 
a 


13.86023 11. 02°3 § 


J0359+509 FSRQ 1.512 12.7973 10.26122 | J0920+446 FSRQ 2.19 12.2577 10. 51733 


J0403-360 FSRQ 1.417 13.7813 20. 14°37? || J0921+622 FSRQ 1.446 11.0372 3.093 


J0423-013 FSRQ 0.916 13.71729 14.28°3 3 || J0927+390 FSRQ 0.695 11.1873? 2.02" 


J0428-379 BL Lac 1.11 13.98105 22. 23*33,78 || Jo927-205 FSRQ 0.348 11.9177 4. 81°58 


-2.54 


J0433+053 BCU 0.033 12.82: ” 4.04'1% | J0937+501 FSRQ 


Y 
Y 
Y 
N 
N 
Y 
Y 
N 
N 
N 
N 
Y 
Y 
N 

J0242+110 FSRQ N 
Y 
Y 
Y 
Y 
Y 
Y 
N 
Y 
Y 
N 
Y 
Y 
Y 
N 0.276 14.33% 16.75236 
Y 


~~ 2 AKA A Ke A KKK Fd KE AK Fd F KH KH ZAREK F KE KEKE KK KK 


J0449+113 FSRQ 2.153 14.0217 25. 5277.8 |] J0958+655 BL Lac 0.367 14. 7970°8 24.61526 
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续 表 1 


Name Type Fermi z logT,,,/K Ôa Name Type Fermi z logT „/K 0 


(1) 2) G) (4) (5) (6) (1) (2) (3) (4) (5) (6) 
J1018+357 FSRQ N 1.229 12552 8253 ||J1504+104 FSRQ 


< 


1.839 13.19222 13.81865 


J1037-295 FSRQ 0.312 14.492 18.73221, ||J1512-090 FSRQ 0.36 14. 157 13.869 


J1041+061 FSRQ 1.264 12.5577 8.21452 |J1516+002 BL Lac 0.052 10.3977? 0.89 


J1043+241 FSRQ 0.56 13. 83° 16.24143 ||J1517-243 BL Lac 0.048 11.47708 1.545 


J1044+809 FSRQ 1.254 13.43122 16. 15*8%8 | J1522-275 BL Lac 1297 92 0.82 


J1048+717 BCU 1.15 14.82:% ° 40.3433 | J1549+026 FSRQ 0.414 13.87183 14.16109 


J1048-191 FSRQ .595 13.3277 13.2474 ||J1550+054 FSRQ 1.422 12.4778 8.751729 


J1058+015 BCU 1.397 11.38772 3.319 


0 
0.89 13.5 12. 09°? % || J1613+342 FSRQ 
0 


J1104+382 BL Lac .03 7.74532, 0.1870 3 | J1625-254 BCU 0.786 13.62225 13.86111 


J1127-189 FSRQ 1.048 23.5715% 6.02'7% | J1626-298 BCU 0.815 14.21208 21.05121 


J1130-148 FSRQ 1.184 13.68111 19.6212% ||J1635+381 FSRQ 1.814 14.06253 27.9178 8 


J1146+399 FSRQ 1.089 13.051532 10.8734 ||J1637+472 FSRQ 0.735 12.5113 6.955 


J1153+495 FSRQ 334 13.05°1% 7. 38°53, ||J1638+573 FSRQ 0.751 13.457) 38 13.92525 


J1159+292 FSRQ 729 13.6775, 13. 29°33 || J1640+397 FSRQ 1.66 13.18 15.9732 


J1159-224 BCU 565 13.2577 11. 05°S 8 || J1642+398 FSRQ 0.593 14.4 21. 6°23 


J1203+480 FSRQ 817 10.2573 1.183 || J1642+689 FSRQ 0.751 12.9471 11.9329 


J1222+042 FSRQ 966 14. 71°72 34. 11°97 | J1658+076 FSRQ 0.621 12.173 5.785 


J1224+213 FSRQ 12.187) 3.465% || J1719+177 BL Lac 0.137 10. 84°37; 1.75749 


158 13.255% 5.96713) || J1727+455 FSRQ 0.717 13.43% 14.3453] 


J1239+075 BCU 4 9.95779 1.153 | J1728+044 FSRQ 0.293 13. 267578 8. 7973 3h 


J1246-257 FSRQ 635 13.967; 4 16. 1373S) || J1733-130 FSRQ 0.902 13.02" 9.481235 


= J1254+116 FSRQ 872 12.4412% 6.6°% | J1743-038 FSRQ 1.057 12.98%) 52 9.961131 


5 © co © co oo os O 5 © 
+. 
Ww 
+ 


J1256-057 FSRQ 536 13.3975 % 8.85 |J1751+096 BL Lac 0.322 13.627) 11.58115 


J1258-223 FSRQ 1.303 12.277? 10.6416% || J1800+784 BL Lac 0. 68 e 9. 7S 
J1305-105 FSRQ 0.286 12.1671 4.47772 | J1806+698 BL Lac 0.046 11.88% 2.1787 
J1310+323 BCU 0.997 13.7812 17. 63°23 | J1824+568 BL Lac 0.663 12.6773} 7.425 


J1316-336 FSRQ 1.21 13.57%) 3! 17. 867°. || J1830+063 FSRQ 0.745 12.397 5.981134 


J1325-430 BCU 0.002 8.08718 0.132 || J1842+681 FSRQ 0.472 11.55 3. 68°73 


J1327+221 FSRQ 1.398 13.58113 20.8778 || J1849+670 FSRQ 0.657 13.43% 12.3639 


J1337-129 FSRQ 539 13.697) 8 13. 29S || J1902+319 FSRQ 0.635 10.82 1. 817933 


J1357+193 FSRQ 72 13.96713, 7.997% || J1911-201 FSRQ 1.119 13.86%) 38 19.46113 


J1415+133 BCU .247 13108 5. 6971) || J1923-210 FSRQ 0.874 13.33 12. 11°88 


J1419+543 BL Lac .153 12.097 2.64 | 1924-292 FSRQ 0.353 14. 17°28 15. 23°38 


J1427-421 FSRQ 1.522 13.08; # 14.23°8% || J1927+739 FSRQ 0.302 12.97" 6.173 


~ 2A + KK K 222 K Z Awe KK Fd AK KK Ke Z AK F AK KEKE KK AK KK Z K 


N 
N 
Y 
Y 
Y 
N 
Y 
Y 
Y 
N 
Y 
Y 
Y 
Y 
N 
N 
Y 
J1229+020 FSRQ Y 
N 
Y 
N 
Y 
Y 
N 
Y 
Y 
Y 
N 
Y 
N 
N 
Y 
Y 
N 


J1446+173 FSRQ 1.024 9.823 0.99% | J2000-178 FSRQ 0.652 14.65112 27. 857252 
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续 表 1 
Name Type Fermi z logT „/K 0 Name Type Fermi z logT,,. /K 6 
(1) (2) (3) (4) (5) (6) (1) (2) (3) (4) (5) (6) 
J2007+404 FSRQ N 1.736 9835 1.19235% |J2225-049 FSRQ Y 140 113 3.063 
J2011-157 FSRQ N 1.178 13.58'F8 21.87°5% | J2229-085 FSRQ Y 156 148109 44.88152 
J2015+371 BCU Y 0.859 13.7771 9% 13.26 | J2232+117 FSRQ Y 1037 14.79*F3 46.865% 
J2025+337 BCU N 0.219 12543 4313 || J2236+284 FSRQ Y 0.79 1455 18. 9177 2 
J2038+513 FSRQ N 1.686 11726 6.46’;% |[J2246-121 FSRQ N 0.632 105322 2. 62°12 
J2123+055 FSRQ Y 1941 13,595 33. 52°17 i || J2258-279 FSRQ Y 0926 12.81553 7. 8710 
J2148+069 FSRQ N 0.999 134163 11.61222 || J2327+096 FSRQ Y 1843 14.46127 42.517} 45 
J2158-150 FSRQ Y 0.672 14.177} 21.19% 3 || J2333-237 BCU N 0.048 11.99% 26310 
J2202+422 BLLac Y 0.069 13.17532 4.152722 |[J2334+076 FSRQ Y 0.401 1253706 7.357393 
J2203+174 FSRQ Y 1.076 10.75°7* 2.2833) | J2348-165 FSRQ Y 0.576 13.6973) 13. 32°33 
J2217+243 BL Lac Y 10.72: 1.85’59 ||J2358-103 FSRQ N 1638 14.38156 34.89: 
J2218-035 FSRQ N 090 18872 58310 
! Fermi detected T 
0.30 Fermi non-detected i 
1 
0.25 
> = 0.20 
3 3 0.15 
= ! = 
0.10 
: 0.05 
0.00 
log Tvar / K 
图 4 亮 温度 分 布 
Fig.4 Brightness temperature distributions 
0.40 
0.35 Fermi non-detected 
0.30 
= 0.25 š 
5 0.20 i 
= 一 


图 5 光 变 多 普 勒 因子 分 布 
Fig.5 Variability Doppler factor distributions 
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利用 Wilcoxon 秩 和 检验 ， 我 们 分 别 比 较 了 R2 焰 变 体 子 类 的 亮 温度 和 多 普 勒 因子 中 位 数 及 误差 
蝎 虎 座 BL 型 天 体 和 平 谱 射 电 类 星体 、 费 米 焰 谈 Table 2 Median brightness temperatures and Doppler 


体 和 非 费 米 炊 变 体 的 亮 温 度 和 多 普 勒 因子 分 布 。 factors of different blazar subclasses and 

每 个 子 样 本 的 亮 温 度 和 多 普 勒 因子 中 位 数 及 误差 a 

见 表 2。Wilcoxon 秩 和 检验 结果 表明 : 2 idles nace Median brightness Median 

(1) 蝎 虎 座 BL 型 天 体 和 平 谱 射 电 类 星体 亮 i temperature/K Doppler factor 

温度 分 布 不 存在 显著 差异 (p=0.12), 但 两 者 的 BL Lacs Bi 4, 1511356 

多 普 勒 因子 分 布 存 芷 显著 差异 (= 0. 006), 平 谱 FSRQs K H W A 11. 247897 

射电 类 星体 平均 具有 更 大 的 多 普 勒 因子 ( 11. a +1.29 +15. 74 
i BCUs 12.92115 5.311152 

4.15: ja 9] 认 为 平 谱 射 电 > 中 
j; RAN x ] AAAH S AURA B Fermi blazars 13. gia 13. 2. 


座 BL 型 天 体 的 观测 视角 (6) 可 能 没有 显著 差异 ， 
但 前 者 具有 相对 更 快速 运动 的 喷 流 ， 至 少 能 解释 NonFemibeas — 1224455 GATS 
这 两 种 次 变 体 子 类 之 间 的 多 普 勒 因子 分 布 差异 。 

(2) 费 米 泡 变 体 和 非 费 米 泡 变 体 的 亮 温 度 和 多 普 勒 因子 分 布 均 存 在 显著 差异 (p 值 分 别 为 1. 73x 
105 和 2. 3x10“)。 相 对 非 费 米 粹 变 体 ， 费 米 炊 变 体 平均 具有 更 大 的 亮 温 度 (13.34? 7 vs 12.44'+4K) 
和 多 普 勒 因子 (13.27*%%, vs 6.17%?)。 这 个 统计 结果 暗示 光 变 体 是 否 能 被 探测 到 伽 马 射线 辐射 可 能 与 
EEDA RO 。 

X JEA IAS I EE TOR ATE 15 GHz 波段 的 光 变 多 普 勒 因子 。 将 文中 样本 和 他 们 的 样 
本 进行 交叉 证 认 ， 发现 具 有 123 个 公共 源 ， 其 中 包括 16 个 蝎 虎 座 BL 型 天 体 ，95 个 平 谱 射 电 类 星体 ， 
12 个 未 知 类 型 的 次 变 体 。 图 6 和 图 7 给 出 了 123 个 公共 源 1 mm 和 15 GHz 的 亮 温 度 和 多 普 勒 因子 比较 
情况 。 在 15 GHz 波段 ， 亮 温度 分 布 的 中 位 数 为 13.97219. 在 1 mm 波段 ， 亮 温度 分 布 的 中 位 数 为 
13.0922, WAKT. Wilcoxon 秩 和 检验 结果 也 表明 ，1 mm 波段 亮 温度 小 于 15 GHz 波段 (p 值 为 
9.47x10”)。 在 15 GHz 波段 ， 多 普 勒 因子 分 布 的 中 位 数 为 12.23152; 在 1 mm 波段， 多 普 勒 因子 分 
布 的 中 位 数 为 9.773442 ， 前 者 大 于 后 者 。Wilecoxon 秩 和 检验 结果 同样 表明 ，1 mm 波段 多 普 勒 因子 平 
均 小 于 15 GHz 波段 (p 值 为 0. 028)。 


logT,,, 15 GHz / K 
Guar_15 GHz 


6 8 10 12 14 16 


logTar_l mm/K ó, 


Var 一 


1 mm 


图 6 1mm 与 15GHz 疲 段 的 亮 温度 比较 。 三 角形 表示 蝎 图 7 1mm 与 15GHz 波段 的 多 普 勒 因子 比较 。 三 角形 表 


虎 BL 型 天 体 ， 圆 点 表示 平 谱 射 电 类 星体 ， 红 色 表 示 蝎 虎 BL 型 天 体 ， 圆 点 表示 平 谱 射电 类 星体 ， 红 

TRAE RRA A, 2k @ 4 F Ht ROM E K ERTI RRA, RERE k KARA 
Fig.6 Comparison of brightness temperatures between 1 mm Fig. 7 Comparison of Doppler factors between 1 mm and 15 
and 15 GHz bands. Triangles represent BL Lac, dots GHz bands. Triangles represent BL Lac, dots represent 
represent FSRQs, red indicates non-Fermi blazars, FSRQs, red indicates non-Fermi blazars, and green 


and green indicates Fermi blazars indicates Fermi blazars 
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文 [16] 和 文 [ 18] 发 现 部 分 耀 变 体 的 射电 亮 温度 可 能 随 观 测 频率 增加 而 降低 ， 多 普 勒 因子 自然 也 
随 频率 增加 而 降低 。 他 们 认为 射电 光 变 多 普 勒 因子 依赖 观测 频率 的 现象 可 能 与 喷 流 加 速 或 喷 流 弯曲 有 
关 。 不 同 频率 的 射电 观测 定位 在 不 同 的 喷 流 区 域 ，1 mm 波段 观测 比 15 GHz 波段 观测 更 靠近 喷 流 内 
区 。 在 喷 流 加 速 模型 下 ， 两 个 区 域 中 喷 流 成 分 运动 速度 的 差异 可 能 导致 mm 波段 和 15 GHz 波段 的 多 
普 勒 因子 不 同 。 另 外 ， 距 流 的 弯曲 现象 导致 在 不 同 射 电波 段 的 观测 视角 不 同 ， 这 也 可 能 导致 1 mm 波 
段 和 15 GHz 波段 的 多 普 勒 因子 不 同 ， 但 文 [16] 和 文 [18] 认 为 这 种 可 能 性 相对 较 低 。 在 本 文中 ， 我 们 
基于 一 个 统计 样本 ,发 现 了 射电 亮 温 度 和 光 变 多 普 勒 因子 对 观测 频率 的 依赖 性 ， 支 持 文 [16] 和 文 
[18 ] 的 部 分 研究 结果 ， 并 且 首 次 将 他 们 的 研究 结果 推 至 亚 毫米 波段 。 但 需要 强调 两 点 : (1 ) 样 本 容量 
较 小 ， 并 且 样 本 由 平 谱 射 电 类 星体 主导 (112/155)。 样 本 中 的 蝎 虎 座 BL 型 天 体 较 少 ( 仅 25 个 ) ， 因 此 
将 来 扩充 样本 (特别 是 加 入 更 多 的 射电 亮 蝎 虎 座 BL 型 天 体 ) 来 验证 文中 结论 很 有 必要 ; (2) 我 们 的 研 
究 是 基于 炊 变 体 子 类 间 的 整体 样本 比较 ， 需 要 注意 到 文 [16] 发 现 射电 光 变 多 普 勒 因子 对 观测 频率 的 
依赖 只 是 对 于 部 分 源 而 言 ， 也 存在 一 些 源 的 多 普 勒 因子 随 着 观测 频率 的 增加 而 增加 ， 或 者 不 明显 依赖 
观测 频率 。 

致谢 感谢 亚 毫米 波 阵 列 的 首席 科学 家 Mark Gurwell 教授 授权 我 们 使 用 光 变 曲线 数据 。 亚 毫米 
波 阵 列 是 美国 史 密 森 天 文 台 和 中 国 台湾 天 文 与 天 体 物 理 研 究 所 的 联合 项 目 ， 由 史 密 森 协会 和 中 研 院 次 
助 。 亚 毫米 波 阵 列 位 于 夏威夷 土著 居民 的 文化 圣地 莫 纳 克 亚 (Mauna Kea) ， 被 特许 在 莫 纳 克 亚 山峰 探 
RFE 
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Estimation of the Doppler Factors for Blazars Using a 
Bayesian Hierarchical Model 


Ning Gongming, Huang Weitian, Mao Lisheng” 


(School of Physics and Electronic Information, Yunnan Normal University, Kunming 650500, China, Email; maolisheng@ynnu.edu.cn) 


Abstract: Based on the long-term monitoring data at 1 mm wavelength from the Submillimeter Array 


(SMA), we used a Bayesian hierarchical model to fit the light curves of 155 blazars and estimated their radio 


brightness temperatures and variability Doppler factors. The brightness temperature and Doppler factor 


distributions between BL Lacertae objects (BL Lacs) and flat spectrum radio quasars (FSRQs) , as well as 


between Fermi and non-Fermi blazars, were compared using the Wilcoxon rank-sum test. The results show that 


(1) FSRQs have higher average Doppler factor than BL Lacs; (2) Fermi blazars have higher average 


brightness temperature and Doppler factor than non-Fermi blazars. After comparing the brightness temperatures 


and variability Doppler factors at 1 mm and 15 GHz wavelengths, we found a tendency of decreasing towards 


higher frequencies for these two parameters. 


Key words: Blazar; Doppler factor; brightness temperature ; Bayesian hierarchical model 


